Avalanche Behavior in an Absorbing State Oslo Model 
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Self-organized criticality can be translated into the language of absorbing state phase transitions. 
Most models for which this analogy is established have been investigated for their absorbing state 
characteristics. In this article, we transform the self-organized critical Oslo model into an absorbing 
state Oslo model and analyze the avalanche behavior. We find that the resulting gap exponent, D, 
is consistent with its value in the self-organized critical model. For the avalanche size exponent, r, 
an analysis of the effect of the external drive and the boundary conditions is required. 
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Self-organized criticality (SOC) refers to the tendency 
of non-equilibrium dissipative systems, with many de- 
grees of freedom driven at a slow rate, to display scale 
invariance without the fine-tuning of any control param- 
eter 01 Q. SOC is present in many open systems where 
the activity, that is, the presence of avalanches, trans- 
ports slope particles through the system to the bound- 
aries where they are eventually dissipated. When there 
is no activity, the system is in a so-called absorbing state. 
The separation of rclaxational and driving timescales is 
achieved by adding a slope unit only when the system is 
in an absorbing state. Thus, the dynamics of the model 
implicitly tune the slope density to values that are asso- 
ciated with the transition between absorbing (inactive) 
and active states. 

In a closed system, an absorbing state (AS) phase tran- 
sition refers to the transition from an absorbing (inactive) 
state to an active state of the system at a critical value 
of a control parameter such as the slope density . 

In retrospect, Tang and Bak's 1988 description of self- 
organized criticality contains the ideas and features of 
absorbing state phase transitions Q . This link was later 
clarified by Vespignani and Zapperi [f| and a recipe for 
transforming AS models into equivalent SOC models was 
devised 6] : Start with a system displaying an AS phase 
transition. When the control parameter is above its crit- 
ical value, there is activity in the system. This activ- 
ity should be coupled via the dynamics to a decrease in 
the control parameter, as for example when the activity 
reaches the boundary. Conversely, when the control pa- 
rameter is below its critical value, there is no activity and 
the system is in an absorbing state. A process, such as 
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the external drive, should increase the control parameter 
by a small amount which may force the system into an 
active state. However, this picture ignores, for example, 
the problem of defining observables common to both AS 
and SOC models as well as the role of finite-size effects. 

In the following, we transform the SOC Oslo model 
with one open boundary into an AS Oslo model with pe- 
riodic boundary conditions. The difficulties in this pro- 
cedure are shortly discussed. Extensive numerical simu- 
lations are analyzed with respect to scale invariance. The 
results lead to a discussion on the effect of the external 
drive and the boundary conditions. 

The model was inspired by an experiment, con- 
ducted in Oslo, on slowly driven rice-piles displaying 
self-organized criticality A one-dimensional lattice 
of length L is characterized by a slope variable, Zi, and 
a critical slope, zf, assigned to each site i = 1,...,L. 
After initialization with z% = and zf drawn randomly 
with equal probability from {1, 2}, the model is updated 
as follows: Driving: A slope unit is added to the leftmost 
site i = 1, such that z\ — » z\ + 1. Toppling: If zi > zf at 
a site i, one slope unit is moved to each of the two near- 
est neighbors, that is, Zi — > Zi — 2 and zi±\ — > Zi±\ + 1 
except when site i = 1 topples, where z\ — > z\ — 2 and 
Z2 — > Z2 + I or when site i = L topples, where zl — > zt, — 1 
and zl-i — > zl—i + 1- After each toppling a new value 
for z\ is chosen randomly with equal probability from 
{1,2}. The activity stops when z% < zf everywhere. The 
model is then driven again. 

After a transient, the slope density, ( = (l/L) Y^,t=i z ii 
fluctuates about a constant value. The avalanche size, s, 
is the total number of topplings after the addition of a 
slope unit. The avalanche size probability density func- 
tion, P(s; L), in a system of size L follows simple scaling 
above a lower cutoff st, 

P(s; L) = as- T G (^) for a > s e , (1) 

where D is the gap exponent (avalanche dimension) and 
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r the avalanche size exponent. The two constants a and 
b are metric factors. The scaling function Q(x) falls off 
sufficiently fast such that all moments of the avalanche 
size probability density function exist in a finite system. 
For n > t—1, the leading order of the nth moment is [9j 

/>oo 

(s n ) (L) = / s n P(s; L)ds = a(bL D ) n+1 - T g n + - ■ • , (2) 
Jo 

where g n depends only on the scaling function Q ■ The 
sub-lea ding terms represent, for example, corrections to 
scaling and the presence of a lower cutoff. 

The numerical estimates of the exponents for the one- 
dimensional SOC Oslo model are very well established 
with D = 2.25(2) and r = 1.555(2) @, ES ELjnde- 
pendent of the exact boundary condition at i = L |18| . In 
the following, we address the key questions: Is it possible 
to recover these exponents, characterizing the avalanche 
behavior and therefore the SOC aspect of the model, in 
its AS counterpart? Is the critical slope density identical 
in both models? 

The recipe mentioned in the introduction can be "in- 
verted" to transform the SOC Oslo model into an AS 
version by imposing periodic boundary conditions such 
that sites i = 1 and i — L are nearest neighbors, that 
is, for all sites toppling (including i — 1 and i = L), 
one slope unit is moved to each of the two nearest neigh- 
bors. The resulting model is translationally invariant, so 
that the original interpretation of the slopes as height 
differences between columns of rice breaks down plQ. 
The slope density increases in steps of 1 /L when the sys- 
tem is driven. If the external drive triggers activity, the 
avalanche will propagate until the systems falls into an 
absorbing state again. Contrary to the SOC Oslo model, 
there is no dissipation mechanism coupled to the activity 
to decrease the slope density in the AS Oslo model. 

Every finite system will eventually fall into an absorb- 
ing state, provided that £ < 2. However, numerically it 
becomes clear that for system sizes L ^> I there exists 
a critical slope density, ( c < 2, above which the absorb- 
ing states become practically inaccessible. Defining the 
activity as the density of sites where Zi > zf, the ac- 
tivity picks up sharply at the critical slope density, £ c . 
However, one fundamental problem with a definition of 
an instantaneous activity is the Abelian nature of the 
Oslo model [l4|, meaning that the order in which sites 
are relaxed is irrelevant. Without an a priori order of re- 
laxation, there can be no a priori microscopic timescale, 
and the temporal behavior of the instantaneous activity 
depends on the choice of microscopic timescale. It will 
be argued below that ( c and all relevant observables can 
be be obtained while avoiding these ambiguities. 

Starting from an empty configuration, £ = 0, slope 
units are added at site i — 1 although, of course, all sites 
are equivalent because of translational invariance. Just 
as in the SOC Oslo model, slope units are added only 
when activity has ceased. Thus, one obtains avalanche 
sizes for slope densities increasing in steps of 1/L. At 



small slope densities, only small avalanches occur and 
until avalanches wrap around the system, the system size 
cannot have any effect on the dynamics. In this regime, 
we find that the behavior of the model essentially resem- 
bles that of the one dimensional BTW model [lj where 
the avalanche size scales like the square of the number of 
slope units added, (CL) 2 . 

At sufficiently large slope densities, the avalanche sizes 
start to deviate from the above behavior. According to 
the original arguments @, one expects that at the AS 
critical point, £ c , the avalanche size probability density 
function follows simple scaling, Eq. Q). Let (s") (£; L) 
denote the nth moment of the avalanche size probability 
density function in the AS Oslo model of size L at slope 
density £. Ignoring corrections to scaling and the exis- 
tence of a lower cutoff in Eq. , one expects that at ( c 
there exist exponents 

7 „ = D(n +1-t) (3) 

such that 

{sn) ^ L) = ab n+1 -'g n forL»l. (4) 

Note that the right hand side is independent of L but 
depends on n, the order of the moment. However, 
away from £ c the ratio on the left hand side depends 
on L. Therefore, for all n and three distinct system sizes 
L\, L2, L3, there exist a unique slope density, £ c , and ex- 
ponent, 7„, such that 

(s")(C e ;£i) _ (s")(C c ;L 2 ) _ (s n )(t c ;L 3 ) 

Tin. 7-7ii r7r> ' \ ) 

L, x 1j 2 L 3 

Graphically, £ c and 7„ are determined by plot- 
ting, for each system size, the rescaled nth moment, 
(s n ) (C; L)/L 7 ", versus the slope density ( and adjust- 
ing the exponent j n until the graphs intersect at a single 
point at £c, see Fig.^ The slope of the graphs of the re- 
sulting exponents, 7„, versus the order of the moment, n, 
determines the gap exponent, D, while the intersection 
with the n-axis gives t — 1, see Eq. J3J). 

Since the moments are measured for slope densities 
increasing in steps of 1/L, the location of the crossings 
requires interpolation between the data points. This in- 
troduces arbitrariness, which can, however, be reduced 
strongly by choosing system sizes which are commen- 
surable with an estimated critical slope density such 
that C C L e N. In a preliminary simulation with L = 
1024,2048,4096,8192 we estimated ( c = 1.73265 which 
suggests L = 1268, 2536, 5072, 10144 as suitable system 
sizes. The numerical estimates presented below refer to 
an exponential interpolation scheme but they are very 
similar to those obtained from linear interpolation. All 
systems were initialized with ( = and gradually filled 
to slightly above ( c , where the average avalanche size be- 
comes extremely large. In each simulation, a configura- 
tion slightly below £ c was recorded and the model driven 
from this configuration 100 times to improve statistics. 
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For every system size, we average over approximately 
2000 independent realizations starting from £ = 0. Thus 
there is a total of approximately 200, 000 correlated re- 
alizations in the neighborhood of £ c . The error bars re- 
ported below are based only on the 2000 independent 
realizations, thereby grossly overestimating the errors. 

The numerical estimate for the critical slope density, 
Co is determined as the slope density at which the graphs 
for different system sizes cross for the rescaled nth mo- 
ment, see Fig. 
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FIG. 1: The logarithm of the rescaled moments 
(a™) {C,;L)/L ln versus the slope density C for n = 1, 2 (open 
circles). The exponential interpolations between the data 
points for increasing system sizes are marked with lines of 
increasing dash length. For the triplet L = 2536, 5072, 10144, 
the lower set of graphs for n = 1 intersect in a single point at 
Cc = 1.732608 with 71 = 2.064 and the upper set of graphs 
for n = 2 intersect at Cc = 1.732609 with 72 = 4.342. 

In principle, each triplet of system sizes produces an 
estimate of Cc for each moment. For example, from the 
crossings of the triplet L = 2536, 5072, 10144 we find 
Cc = 1.73261(5) for n = 1 and Cc = 1.73261(5) for 
n = 2, see Fig. □ For the triplet L = 1268, 2536, 5072 
the resulting estimates are Cc = 1.73262(8) for n = 1 and 
Cc = 1.73259(7) for n — 2. Ignoring errors, all triplets 
and moments yield 

1.73257 < Cc < 1.73262 . (6) 

There is no established systematic way to determine the 
value of Cc in the limit L — > 00 but from the small change 
in the estimate as the system sizes are increased, it seems 
reasonable to estimate the critical slope density by av- 
eraging over all moments n = 1,2,3,4 and all triplets 
yielding Cc = 1.73260(2). 

Since crossings rely on an appropriate choice of 7„ , es- 
timates of the critical slope density Cc and "f n go hand 
in hand, see Fig. ^ Plotting 7„ versus the order of 
the moment, n, produces according to Eq. @ an esti- 
mate for D and r. This procedure can be performed for 
every triplet of system sizes. For L = 1268, 2536, 5072 



we find D = 2.22(23) and r = 1.07(18) while for 
L = 2536, 5072, 10144 the resulting exponents are D — 
2.28(23) and r = 1.09(13), see Fig.EJ 




FIG. 2: The estimated exponents 7„ versus the order, n, 
of the moment for the triplets L — 1268, 2536, 5072 (open 
squares) and L = 2536,5072,10144 (open circles). Assum- 
ing simple finite-size scaling it follows from Eq. that the 
slope of the graphs is the gap exponent, D, while the in- 
tersection with the n-axis is r — 1, see inset (dotted lines 
indicate magnified region). Linear regression yields D — 
2.22(23), t = 1.07(18) for the first triplet (dashed line) and 
D = 2.28(23), r = 1.09(13) for the second triplet (solid line). 

The errors of these exponents were estimated by stan- 
dard error propagation. Just as for the estimate of Cc, the 
procedure only provides exponents for a particular set of 
system sizes. However, even the preliminary simulations 
with L = 1024, 2048, 4096, 8192 are fully compatible with 
the exponents D and r reported above. Based on all four 
possible triplets of system sizes, it therefore seems justi- 
fied to estimate D = 2.25(8) and r = 1.08(5), in the limit 
L — > 00. 

Table [I] summarizes the numerical estimates for the 
critical slope densities and exponents for the AS and 
boundary driven SOC Oslo model. 





Quantity 


Oslo model 




D 


T 


7i 


AS 


1.73260(2) 


2.25(8) 


1.08(5) 


2.07(6) 


SOC, boundary driven 


1.7326(3) 


2.25(2) 


1.555(2) 


1 


SOC, bulk driven 


1.734(2) 


2.25(3) 


1.10(3) 


2 



TABLE I: The critical slope density Cc, the gap exponent 
D, the avalanche size exponent r, and the exponent 71 in 
the AS Oslo model (all errors based on averaging over the 
four triplets), SOC Oslo model driven at the boundary 
HI [3 or in the bulk For SOC Oslo models, a simple 

conservation argument in the stationary regime determines 

We are not aware of a systematic study of the aver- 
age slope density (C) (L) in the SOC Oslo model. How- 
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ever, results from simulations of systems sizes L = 
1024, 2048, 4096, 81 92 are consistent with ( c - (() (L) oc 
L~ x with x ~ 0.7 and Cc = 1.7326(3), where the error 
bar is based on visual inspection. Similarly for the bulk 
driven model (see below) we find £ c = 1.734(2). 

It is surprising how well the non-universal slope den- 
sity of the SOC Oslo model is reproduced by the AS 
model. This, however, is exactly what is predicted by 
the simple mechanism put forward by Vespignani and 
Zapperi 

The avalanche size exponent t = 1.08(5) in the AS 
Oslo model is inconsistent with the value r = 1.555(2) 
in the SOC Oslo model. However, the gap exponent 
D = 2.25(8) in the AS Oslo model is consistent with the 
value D = 2.25(2) reported in the literature for the SOC 
Oslo model [3, [HA 113 • While the former result ques- 
tions the one-to-one correspondence between the SOC 
and AS Oslo model, one might understand this discrep- 
ancy as follows: The exponent actually characterizing 
the Oslo model is the gap exponent D. The gap expo- 
nent is deeply rooted in the model and is also present in 
the corresponding field theory |l2t Il6|. In contrast, in 
the SOC Oslo model, the avalanche size exponent r is 
determined by the gap exponent D via the scaling rela- 
tion Eq. J3J), since 71 can be derived by simple conser- 
vation arguments in the stationary regime; driving the 
SOC Oslo model at any fixed, absolute position, such as 
i = 1 in the boundary driven SOC Oslo model described 
above, leads to 71 = 1, implying D(2 — r) = 1. Driving 
the model randomly in the bulk with slope units or at 
any relative position, such as i = L/2, leads to 71 =2, 
implying D{2 — r) = 2. The latter scaling law produces 
t = 1.111(8) from D = 2.25(2), well compatible with the 
result found in the AS Oslo model. In fact, the expo- 
nents D = 2.25(3) and r = 1.10(3) have been reported in 
the literature for a variant of the bulk-driven Oslo model 
|l5j| . but interpreted as evidence against the bulk-driven 
SOC Oslo model being in the same universality class as 
the boundary driven SOC Oslo model. 

Contrary to the SOC Oslo models, there exists no 



conservation argument to determine the value of 71 in 
the AS Oslo model. However, because of the periodic 
boundaries, the AS Oslo model seems to be related to 
the bulk-driven SOC Oslo model rather than the orig- 
inal boundary-driven SOC Oslo model. Therefore, the 
discrepancy in the avalanche size exponent, t, in the AS 
Oslo model with respect to the SOC Oslo model driven 
at any fixed position might very well be caused by the 
difference in the external drive and the resulting scaling 
relations. In fact, the AS Oslo model analyzed above 
seems to be a perfect AS version of an SOC Oslo model 
driven in the bulk at, say, i = L/2. 

Since the SOC Oslo model is equivalent to an inter- 
face depinning model 0, , the latter should also have 
an equivalent AS version, see discussion in Ref. [13 • It 
is also conjectured that the train model is in the same 
universality class as the Oslo model [l2j. Whether im- 
posing periodic boundary conditions on the train model 
will transform it into its AS analogue is an intriguing 
question which is beyond the scope of this brief report. 

In conclusion, we have transformed the SOC Oslo 
model into an AS Oslo model. We have numerically 
determined the critical slope densities and obtained nu- 
merical estimates for the gap exponent, D, and avalanche 
size exponent, r, characterizing the avalanche behavior 
assuming simple finite-size scaling. The critical slope 
densities and the gap exponents are identical in the two 
models. However, the avalanche size exponents are differ- 
ent. The question remains whether one can construct a 
"proper" AS Oslo model, corresponding to the boundary 
driven SOC Oslo model. 
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